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NONLINEAR FLEXURAL VIBRATIONS OF 
THIN- WALLED CIRCULAR CYLINDERS 


By David A. Evensen 
Langley Research Center 

SUMMARY 

The nonlinear flexural vibrations of thin-walled circular cylinders are analyzed by 
assuming two vibration modes and applying Galerkin's procedure. This procedure 
results in two coupled nonlinear differential equations for the modal amplitudes. Only 
one of these differential equations contains a forcing function, since the applied loading 
is so chosen that only one mode (the driven mode) receives external excitation. 

Approximate solutions for the modal response are obtained by the method of aver- 
aging. One such solution involves the response of only the driven mode; a subsequent 
stability analysis shows that this single-mode response is unstable for certain combina- 
tions of amplitude and frequency. In this unstable region, it is necessary to consider the 
coupled-mode response involving both modes, and an approximate solution is presented 
for this problem. 

The analysis is in qualitative agreement with the available experiments on the non- 
linear vibrations of cylindrical shells, and it exhibits several features that are thought 
to be characteristic of nonlinear flexural vibrations of thin shells of revolution. 

INTRODUCTION 

Current design of missiles and launch vehicles relies extensively on the use of 
thin-walled cylindrical shells as the primary structure. During powered flight, these 
cylindrical structures are often caused to vibrate to large amplitudes. This problem of 
large- amplitude vibrations has given rise to a number of theoretical studies of the non- 
linear vibration of thin cylindrical shells. (See refs. 1 to 4, for example.) The basic 
approach used in these studies is to assume the shape of the deflection in space, that is, 
the shape of the vibration mode, and then to derive a corresponding amplitude-frequency 
relationship. 

In a pioneering work, Reissner (ref. 1) isolated a single half-wave or lobe of the 
vibration mode and analyzed it. His results indicated that the nonlinearity could be 
either of the hardening (frequency increasing with amplitude) or softening (frequency 



decreasing with amplitude) type, depending upon the geometry of the lobe. Reissner's 
work was followed by a similar study by Chu (ref. 2). 

Chu employed the same assumed mode shape as reference 1 but analyzed the 
problem somewhat differently. His analysis indicated that the nonlinearity was always 
of the hardening type, in which case the frequency increases with the amplitude of 
vibration. 

One reason for the different results of references 1 and 2 was discovered by 
Cummings (ref. 3). Cummings employed a Galerkin procedure and found that the results 
varied with the region of integration. When the integration was carried out over a single 
half-wave or lobe, the results obtained were the same as those of Reissner. When the 
integration was extended over the entire shell, the results were similar to those of Chu 
(ref. 2). Thus, it appeared that Reissner's results were characteristic of curved panels, 
whereas Chu's calculations were apparently applicable to complete cylindrical shells. 
Chu's calculations subsequently received further confirmation from a study by Nowinski 
(ref. 4). 

A striking feature of the results from references 2 to 4 is the conclusion that the 
nonlinearity always causes an increase in frequency with amplitude. On the other hand, 
recent experiments (refs. 5 to 7) indicate that the vibration frequency generally decreases 
with amplitude. This difference between the trends indicated by theory and experiment 
led to a reexamination of the previous analyses and thus originated the present study. 

A review of the previous analyses revealed that references 1 to 3 do not satisfy an 
important continuity restraint on the in-plane circumferential displacement. Conversely, 
reference 4 satisfies the circumferential restraint but does not maintain zero transverse 
deflection at the ends of the cylinder. In order to satisfy the continuity condition and 
still have no deflection at the ends of the cylinder, the present report employs a deflec- 
tion shape which is a modification of that used in references 1 to 3. Unfortunately, this 
modified mode shape is not moment-free at the ends of the cylinder. Thus, the mode 
shape used herein has boundary conditions that lie somewhere between simply supported 
and clamped ends. ! 

For the most part, references 1 to 4 have concentrated on the problem of free 
vibrations of a single vibration mode. Cummings (ref. 3) formulated a two-mode forced- 
vibration analysis but gave few quantitative results and was primarily concerned with 
the weak coupling problem. For weak coupling pf the modes, each mode can vibrate 
independently. On the other hand, the present analysis treats a case of strong modal 
coupling, in which both modes vibrate with comparable amplitudes. The two modes 
involved are classified as a driven mode and its companion mode, both of which have the 
same natural frequency. Related studies on nonlinear vibrations of circular rings 
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(ref. 8) and circular plates (ref. 9) have shown the need for including these coupled 
modes in the study of nonlinear forced vibrations of axisymmetric structures. 

In summary, then, the present study involves both free and forced vibrations of a 
thin-walled circular cylinder. For forced vibrations, coupled motions involving a driven 
mode and its companion mode are analyzed. The assumed mode shapes satisfy the cir- 
cumferential continuity condition, but the longitudinal boundary conditions are somewhere 
between simply supported and clamped. In-plane inertia effects are neglected, and 
elastic vibrations are assumed. The nonlinearities considered are of the geometric type, 
originating in the strain-displacement relations. Many features of the analysis are sim- 
ilar to a previous study on rings (ref. 8) and are thought to be characteristic of nonlinear 
vibrations of thin shells of revolution. 
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a mn amplitude coefficient in equation (4a) 


A mn (t)»B mn (t) coefficients of the assumed deflection shape 
A(t),B(t) slowly varying amplitudes 

A,B average values (over one period) of A(r) and B(t), respectively 

C(t) generalized force; as defined in equation (9a) # 

D bending stiffness, — ^ 

12(1 - v 2 ) 

E Young's modulus 
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F(x,y,t) 


stress function 
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stress-function amplitude in equation (4b) 
nondimensional amplitude of the applied loading 
cylinder wall thickness 
Mathieu equation coefficients 
length of the cylinder 

number of axial and circumferential waves, respectively 
in-plane stress resultants 
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matrices used in the stability analysis; see equation (B17) 


radial loading applied to the surface of the cylinder 


spatial distribution of the radial load; see equation (11) 


radius of the cylinder 


generalized force; see equation (9b) 


time 


axial, circumferential, and radial displacements, respectively, of a point 
on the median surface 


stability variable, defined by U(0) = 0 exp 
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x, y axial and circumferential coordinates on the median surface, respectively 

z radial coordinate 

a = mir/L 
/3 = n/R 

y, 6,e nonlinearity parameters defined in equations (13) and (14) 

€ xy> e x> e y median surface shear, axial, and circumferential strains, respectively 

£ c (t),? s (t) nondimensional generalized coordinates associated with A mn and B mn , 
respectively 


9 = Qt 




v Poisson's ratio 


£ aspect ratio, 

p mass density 


t nondimensional time, co mn t 


4>,ri small perturbations in the amplitudes of the driven mode and the companion 

mode, respectively 

4>i,(p2> r li> r l2 stability variables 

(P,<Pq column vectors involving <f > 2 , and tj 2 ; see equation (B17) 

co vibration frequency 

<x>l experimental linear vibration frequency 

w mn calculated linear vibration frequency, defined by 

0,2 = E f I 4 e(| 2 + l) g 

mn + i)2 I2(l - i/2) 
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nondimensional frequency, 


a? 

w mn 


q 

V 4 biharmonic operator, V 4 = -^-rr + 2 — | ~ 

9x 4 dx 2 dy z 9y 4 

A bar over a quantity indicates that it is an average value, where the average is taken 
over one period of the motion. 


THEORY 

The governing equations are first presented and briefly discussed. The problem 
is reduced to one involving ordinary differential equations by applying Galerkin's method. 
Approximate solutions for the nonlinear vibration of both a driven mode and its compan- 
ion mode are obtained by the method of averaging. Stability of the response is also 
discussed. 


Governing Equations 

The well-known approximations of Donnell’s shallow- shell theory for thin-walled 
circular cylinders as exemplified in reference 2 result in three equations of motion. 
These equations can be combined to give 


TXT7 4„ . a 2 W „ . 1 9 2 F . [ 9 2 F 9 2 w o 9 2 F 9 2 w . 9 2 F 9 2 w\ 

Dv w + P h ^ =q + __ + ^_ J _ r .2 S8y5 r^ + _ 2 -_ ;f j 


and 


_Lv 4 F= 1 g2w , ( 9 2 W 9 2 W 

Eh ’Rax 2 \9x 9y/ ax 2 ay2 


( 2 ) 


The coordinate system and shell geometry 
are shown in figure 1. It should be noted that 
equations (1) and (2) are limited by various 
approximations made in their derivation. For 
example, they do not contain the effects of in- 
plane inertia; consequently, they are limited to 
motions that are primarily flexural. In addition, 
only the predominant nonlinear terms have been 
included in deriving these equations. The non- 
linearities in the bending terms have been 
neglected, whereas the nonlinearities in 



Figure 1.- Shell geometry and coordinate system. 
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stretching have been estimated by using the following midsurface strain-displacement 
relations: 


e 


x ~ 


8u l/ &wY^ 
ax 2\dxJ 


_ av _ w 1 / aw\ 
? dy R 2\8y ) 


_ 8u ay aw aw 

xy 8y 8x 9x 9y 

Donnell’s approximation « lj has been used, and this assumption limits the analysis 

to circumferential wave numbers greater than about 5. The usual thin- shell assumption, 
2 

« 1, has been employed in the derivation, and transverse- shear and rotary-inertia 
effects have been neglected. 

The boundary conditions appropriate to equations (1) and (2) vary with the problem 
to be analyzed. However, for a complete cylindrical shell, such as is considered in the 
present problem, it is apparent that the displacements, slope, moments, shears, and 
stresses must all satisfy continuity conditions of the form 



w(x,y,t) = w(x,y + 2irR,t) 
v(x,y,t) = v(x,y + 27rR,t) 


(Note that it is not sufficient that the stress function F be continuous and periodic in 
the y-variable.) 

The preceding continuity requirements place a restriction on the choice of mode 
shapes which can be used in approximate analyses; the mode shapes are discussed in the 
following section. 


Reduction to Ordinary Differential Equations 

Equations (1) and (2) can be solved, approximately, by choosing a vibration mode 
for w, solving equation (2) to find F, and then using Galerkin's method on equation (1). 
Throughout the calculations, w and F must be examined to verify that the necessary 
continuity requirements (eq. (3)) are met. 

Choice of vibration modes .- For free vibrations of a freely supported cylinder, 
the linearized versions of equations (1) and (2) have the solution (ref. 2): 

w(x,y,t) = a mn cos ^ sin cos a> mn t (4a) 


7 


(4b) 


F(x,y,t) = f mn cos ^ sin cos o> mn t 

For nonlinear vibrations of a freely supported cylinder, a deflection of the form 

w(x,y,t) = A mn (t)cos sin (5) 

was used in references 1 to 3. However, as discussed in reference 5, equation (5) 
together with the corresponding stress function will not satisfy the necessary periodic 
continuity constraint on the circumferential displacement v. This difficulty is counter- 
acted by replacing equation (5) with 

w(x,y,t) = A mn (t)cos sin ^ + ^A mn (t)sin (6) 

where the term in the brackets is included to satisfy the continuity requirement on v 
(eq. (3)). Note that equation (6) gives w = 0 at the ends of the cylinder but the axial 
bending moment is not zero there. 

Finally, previous studies on circular rings (ref. 8) and circular plates (ref. 9) 
have shown the necessity of including a companion mode in the nonlinear analysis. This 
additional mode has the same linear vibration frequency and the same (m,n) numbers as 
do the modes given in equations (4) to (6); hence, the term "companion mode" is used to 
describe it. When the companion mode is included, equation (6) becomes 

2 

w(x,y,t) = A mn (t)cos ^ + B mn (t)sin ^Jsin 2^ + 2_ A mn 2 (t) + B mn 2 (t) sin 2 222 ( 7 ) 

Equation (7) represents the deflection modes assumed in the present report; these mode 
shapes are similar to those which apply to rings (ref. 8) and cylindrical shells of infinite 
length. 

Application of Galerkin's met hod.- Before the Galerkin procedure can be applied, 
the stress function F must first be found. This function is determined by substituting 
the assumed deflection (eq. (7)) into equation (2) and solving for the particular solution. 
This procedure gives 

F(x,y,t) = aj (A mn cos /3y + B mn sin /3y)sin ox 

" a 2 (A mn 2 - B mn 2 )cos 20y - a 3 A mn B mn sin 2/3y 

+ a 4 (A mn 2 + B mn 2 )(A mn cos /3y + B mn sin |3y)sin 3ox (8) 

where a = mir/L, /3 = n/R, and the coefficients a 4 to a 4 are defined in the section 
"Symbols." 
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At this stage in the analysis, w and F were examined to make certain that all 
the necessary continuity requirements were satisfied. (See appendix A for a further 
discussion of this point.) With w and F established as appropriate for a complete 
cylinder, equations (7) and (8) are substituted into equation (1), which is approximately 
satisfied by using Galerkin's procedure. This procedure yields two coupled nonlinear 
differential equations for A mn (t) and B mn (t). In nondimensional form, these equa- 
tions are 



and 


IF 


+ r 

S s 8 Ss 


d 2 f 

• s s 
>s dr 2 






(9b) 


where f? and £ are vibration amplitudes, given by 
c s 



(10a) 


and the dimensionless time r is defined by 


T 


mn 


(10b) 


where o) mn is the linear vibration frequency. 

The functions C(r) and S(r) in equations (9) represent the generalized forces 
on the two modes and are given by integral expressions involving q(x,y,t), &w/3A mn , 
and 3w/3B mn . The external loading q(x,y,t) is assumed to be fixed in space and 
harmonic in time, so that 

q(x,y,t) = Q(x,y)cos o>t (11) 
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The function Q(x,y) is assumed to be symmetric with respect to y and to have zero 
average value. In this case S(r) is identically zero and 

C(r) — Gjjju cos £2 t (12a) 

where J2 = - a ) — is a nondimensional frequency and the coefficient G m n is given by 
w mn 

G mn = \ o S' S' Q(x,y)cos ^ sin 2^2 dy dx (12b) 

(LR 7 r)ph 2 w mn 2 J 0 J 0 R L 


Note that every nonlinear term in equations (9) contains the parameter e , which is 
defined by 


Thus, e is the basic nonlinearity parameter in the problem; when e goes to zero the 
vibrations become linear. The other parameters which influence the nonlinearities are 
y and 6; they are defined by 

t 4C 1 1 £ 1 


(| 2 + if 16 12(l - if 

I 4 , €(4 2 + lP1 
[{2 + l f 12(l -y 2 ] 


1 i * 1 — 

_ 3| 4 |({ 2 + l) (9S 2 tl). 

" 16 [ i 4 

(|2 + l) 2 12(l-^) 


where ^ is the aspect ratio, given by 


t _ nur/L _ Circumferential wavelength 
n/R ~ Axial wavelength 


In the derivation of equations (9), it should be noted that the weighting functions 
8w / 8 A m n an d &w/3B mn have been used in the Galerkin procedure. This choice of 
weighting functions makes the resulting differential equations for A mn and B mn 
equivalent to those obtained by using the Rayleigh- Ritz method on this problem. (See 
ref. 10 in this regard.) 


1 1 ■ II !■ 



Furthermore, note that the specification of C(r) = G mn cos Sir as the forcing 

function leads naturally to classifying the displacement associated with A mn (t) and 

cos ^ sin S122E as the driven mode. The companion mode is associated with B mn (t) 
K Li 

and sin ^ sin ; and since S(t) = 0, the companion mode is not directly excited by 
K Li 

the forcing function. 


Approximate Solutions to the Nonlinear Equations 

Inasmuch as C(r) and S(r) have been specified, equations (9) can be solved 
approximately by the method of averaging (ref. 11). Such solutions will be presented for 
vibrations involving (a) only the driven mode, and (b) both the driven mode and its com- 
panion mode. 

Single-mode response .- Because S(r) = 0, a possible solution to equation (9b) is 


that £_(t) is identically zero. 

S 


In this case, only £ (r) vibrates, and gives rise to the 


single- mode response. Using £„(r) = 0 and ?„(t) = A(t)cos S2t in equations (9) and 

S C 

applying the method of averaging gives the approximate solution 


£ c (t) = A cos 

C S M = o 

where the average amplitude A is computed from 

(l - 0 2 )a- f! S2 2 A 3 - A 3 + | e 2 6A 5 = G 
16 4 8 


mn 


(16) 


(17) 


for given values of e, G mn , and the other quantities, 
as a function of is given in figure 2 for e = 0.01 
comparison, e =0.01 corresponds to n = 10- 
and h/R = 0.001, and £ = 0.1 corresponds 
to m = 1, n = 10, and L/R = it. 

The solution for free vibrations involving 
a single mode can be obtained from equa- 
tion (17) by putting G mn = 0. This solution 
gives the one-mode "backbone" curves shown 
in figures 3 and 4. These curves were com- 
puted from the expression 


A typical response curve of | A| 
and £ = 0.1. For purposes of 


O 2 = 


1-^A 2 




+ ^®A 4 
+ 8 A 


(18) 



and give the amplitude-frequency relation for 
free vibrations of a single mode. 


Figure 2.- 


Single mode response. Forced and free vibrations, 

£ =(t -) 2 = 001 ' * = T7m =0 - L 
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Figure 3.- Influence of large amplitudes on vibration frequency for various values of e. Free vibrations-, one mode; 

£, = 1/2 and 2; v = 0.3. 


The stability of the approximate solution given by equation (16) was examined by 
the usual techniques of perturbation analysis, as described in appendix B. The resulting 
Mathieu-Hill stability equations (ref. 12) indicate that within terms of order e^: 

(a) Perturbations of ? are unstable within the area bounded by (eq. (B8)) 

V 

(b) Perturbations of £ are unstable within the region (eq. (B9)) 

s 

#(H 

(c) Both types of perturbations are unstable in narrow regions near 
S2 = 1/2, 1/3, ... 

The first instability region (a) coincides with the locus of vertical tangents to the 
response curves and indicates the well-known jump phenomena. The narrow areas (c) 
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(ttRAi) 

(L/m) 


Figure 4.- 


influence of large amplitudes on vibration frequency for various aspect ratios. Free vibrations; one mode; 


£ = 



1.0; v = 0.3. 


near £2 = 1/2, 1/3, . . ., denote possible ultraharmonic responses. The remaining 
region, (b), indicates the area in which the companion mode is parametrically unstable 
because of nonlinear coupling with the driven mode. If adequate solutions are to be 
obtained in region (b), it is necessary to consider motions where both modes vibrate. 


Co upled-mode response .- When both £ g (r) and £ c (t) oscillate, the method of 
averaging (ref. 11) gives the approximate solution 




= A cos £2 t 
= B sin £2 t 


(19) 


where A and B satisfy the following equations: 
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(l - S2 2 )a+ a(b 2 - A 2 ) - ^3A 2 + B 2 ) + a(5A 4 + 2A 2 ! 2 + B 4 ) = G, 

16 ^ o 


(20a) 


mn 


(l - S2 2 )b+ b(a 2 - B 2 ) - ^(3B 2 + A 2 ) + ^ b(5B 4 + 2A 2 ! 2 + A 4 ) = 0 


(20b) 


Note that one possible solution to equations (20) is that B = 0. In this case, equa- 
tions (19) and (20a) revert to equations (16) and (17), respectively. Furthermore, inspec- 
tion of equation (20b) reveals that real, nonzero values of B exist only if £2 satisfies 
the condition 

S2 < 1 + S£_(| -V) +o(e 2 A 4 ) 

which agrees with the preceding single- mode stability results, region (b). 

When A and B are not zero, they can be computed by solving equations (20) 
simultaneously. A typical coupled-mode response curve for this case is shown in fig- 
ure 5, for e = 0.01 and £ = 0.1. 



Figure 5.- Coupled-mode response. Forced vibrations; e = 



= 0.01; £ = 



0.1; G mn = o.l; V = 0.3. 
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The case of free coupled- mode vibrations can be obtained by putting G mn = 0 
and A = B in equations (20). This procedure yields 

S2 2 = 1 - eyA 2 + e 2 6A 4 (21) 

which was used to compute the two-mode backbone curves shown as the dashed curves 
o-n in figure 5. 

The stability of the approximate solution given by equations (19) was investigated 
by using small perturbations and applying the method of slowly varying parameters. 
Details of the stability analysis are given in appendix B. Based upon this analysis, 
branches a-b and b-c of the response curves in figure 5 appear to indicate stable 
responses, as does the segment f-g. The other branches of the response curves, namely, 
c-d, e-f, and g-h-i-j in figure 5, appear to be unstable. 

DISCUSSION OF RESULTS 
Single-Mode Results 

A typical single-mode response curve is shown in figure 2, which illustrates a non- 
linearity of the softening type. This behavior is typical for a vibration mode that involves 
low values of e and £. Such cases often occur for a shell that is relatively long and 
thin-walled. For larger values of e , the nonlinearity is stronger, and for some values 
of £ the nonlinearity is of the hardening type. 

The manner in which the parameter e controls the strength of the nonlinearity is 
shown in figure 3. Linear vibrations correspond to e = 0, and increasing e causes the 
vibrations to become more and more nonlinear. This result is apparent from equa- 
tions (9) and (17), which show that e is a multiplying factor in every nonlinear term. 
From a physical standpoint, small values of e correspond to very thin-walled cylinders 
and/or a low number of circumferential waves, whereas large e signifies thick-walled 
cylinders or a high number of circumferential waves. 

Another important parameter is the aspect ratio £, which governs the character 
of the nonlinearity (that is, softening or hardening). The way in which the type of non- 
linearity varies with the aspect ratio is shown in figure 4. Actually, the parameters e, 
y, and 6 in equations (9) and (17) combine to determine the type of nonlinearity present; 
but both y and 6 are strong functions of £. Small values of £ (say £< 1/2) result 
in small values of y and 6, and the corresponding terms in equation (17) are relatively 
insignificant for small vibration amplitudes. In this case, the vibrations exhibit a soft- 
ening nonlinearity for small amplitudes. As the amplitude continues to increase, how- 
ever, the e 2 6A 5 term in equation (17) finally dominates the other terms and causes an 
eventual hardening nonlinearity. 



For large values of £ (say £ > 2) the y and 5 terms in equation (17) domi- 
nate the nonlinearities immediately, and a corresponding hardening type of nonlinearity 
is observed for all amplitudes. The remaining area, where i < £ < 2, represents a 
transition region in which the parameter y experiences a change in sign. Note that 
from its definition, small values of the aspect ratio £ correspond to short circumferen- 
tial and long axial wavelengths, whereas the reverse is true for large values of 


Of particular interest is the limiting case when the aspect ratio tends to zero. 
Such a case occurs when the length L tends to infinity. When free vibrations of a 
single mode are examined as L — °°, both y and 6 vanish and equation (9a) reduces 
to 


d 2 ^ 


dr 


+ £ + — ? 

2 b c 8 s c 


d\ 


fas. 


-r* 


dT 2 


dT, 


= 0 


( 22 ) 


This equation gives rise to a softening type nonlinearity and bears a strong resemblance 
to results obtained for inextensional nonlinear vibrations of rings (ref. 8). 

For comparison purposes, Chu's results (ref. 2) have been plotted in figure 4. The 
theory of reference 2 indicates that the nonlinearity is always of the hardening type; sim- 
ilar results were obtained in references 3 and 4. It should be noted that the results of 
references 2 to 4 indicate a symmetric dependence on the aspect ratio |, whereby the 
curves for £ = • • • coincide with those for £ = 2, 4, 8, . . ., respectively. 

Such a symmetric dependence on the aspect ratio occurs for the nonlinear vibration of 
simply supported flat rectangular plates (ref. 13), as would be expected from the physical 
symmetry of the plate problem. 

On the other hand, the problem of a thin-walled circular cylinder is physically non- 
symmetric, since the cylinder is curved circumferentially but not axially. In the present 
problem then, it appears that the calculations should exhibit a nonsymmetric dependence 
on the aspect ratio £. Such a nonsymmetric dependence on £ is borne out by the 
results of the present analysis, as figure 4 shows. 


Coupled-Mode Results 

A plot of a typical response involving a driven mode and its companion mode is 
shown in figure 5. The results are analogous to those obtained for nonlinear vibrations 
of rings (ref. 8) and for nonlinear vibration absorbers (ref. 14). Along the a-b portion 
of the curves (fig. 5), B is zero and only the driven mode responds. The segments 
b-c-d and e-f-g indicate vibrations in which both modes have comparable amplitudes; 
note that they are adjacent to the two-mode "backbone" curves, o-n. Along g-h-i-j, the 
response of the companion mode is much greater than that of the driven mode; note that 
the response of B as a function of and the one-mode backbone curve (o-h-i-j) in 


16 



figure 5 then coincide. This result is analogous to the vibration- absorber response in 
which the driven mass experiences very little motion while the absorber mass vibrates 
with large displacements. 

Points of vertical tangents to the response curves occur at c,f,g, and h, and the 
stability analysis indicates that the segments c-d, e-f, and g-h-i-j of the response curves 
are unstable. Instead of responding along the segment g-h-i-j, beyond the point g the 
vibrations apparently revert to the one-mode case, with B = 0 and A given by equa- 
tion (17). In some cases, a gap in the solutions may occur, and then both the one-mode 
solution (eq. (16)) and the two-mode solution (eq. (19)) are unstable. Similar unstable 
gaps have been observed in related problems, and analog computer studies indicate that 
nonsteady vibrations with rapidly changing amplitudes occur in these regions. (For 
example, see refs. 8, 15, and 16 in this regard.) 

Comparison With Available Experimental Results 

Experimental results for the nonlinear flexural vibrations of thin-walled circular 
cylinders are rather scarce. Reference 6 gives nonlinear vibration data for cylinders 
both with and without an internal liquid. The results exhibit a softening type of non- 
linearity, and for an empty cylinder the nonlinearity observed was very slight. Similar 
qualitative observations were reported in reference 5, and some recent quantitative 
results have been given by Olson, reference 7. 

Olson's experimental data for forced vibrations are shown by the circles and 
dashed lines in figure 6. The solid line in figure 6 is the free-vibration response curve 



Figure 6 .- Comparison with experimental results, e = = 3.025 x 10"3 ; 

E =^ 7 ^ = 0.1635; v = 0.365. 
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calculated from equation (18) with values of e , y, and 6 that correspond to Olson's 
experiment. Note that the experimental results are nondimensionalized with respect to 
the experimental linear frequency a> L , whereas the analytical results are divided by the 
calculated frequency w mn . Both theory and experiment indicate a nonlinearity of the 
softening type for this case, whereas the results of references 2 to 4 indicate a hard- 
ening type of nonlinearity for all cylinders. 

Olson also detected a "double frequency contraction" at the nodes of cos 

R. 

Because of the presence of the bracketed term in equation (6), the present analysis pre- 
dicts such a double- frequency effect. For sufficiently large amplitudes, nonsteady 
vibrations were also observed (ref. 7); this behavior may have been due to the presence 
of the companion mode. 

Vibrations involving both a driven mode and its companion mode have been found 
for ring vibrations, and traveling-wave responses can result (ref. 8). In the present 
case, the theory yields a traveling wave if the amplitudes A and B are identical. 
When A = B, the deflection becomes 

— = A.(cos Or cos ^ + sin fir sin ^)sin + - — (cos 2 £2t + sin 2 S2r)sin 2 

h \ R R/ L 4R ' 7 L 

= A cosgZ - fir)sin ^ + || A 2 sin 2 (23) 

where the first term represents a circumferentially traveling wave. Traveling- wave 
responses have been observed in the vibration studies of reference 17 and in recent 
flutter tests on thin cylindrical shells (ref. 18) . 

CONCLUDING REMARKS 

The nonlinear flexural vibrations of thin-walled circular cylinders were analyzed 
by choosing vibration modes and applying Galerkin's method. Although only one mode 
was directly driven by the forcing function, it was necessary to include two vibration 
modes in the calculations. Two modes were required because under certain conditions 
nonlinear coupling caused the companion mode to respond and participate in the motion. 
For other conditions, the single-mode response was sufficient. In both cases (single 
and coupled modes) the degree of nonlinearity is dependent upon wall- thickness — radius 
ratio and number of circumferential waves, and the type of nonlinearity is determined by 
the aspect ratio £. Small values of £ generally result in softening characteristics, 
and large values of £ give rise to hardening effects. 
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The single- mode calculations are in qualitative agreement with the available test 
results, which indicate a slight nonlinearity of the softening type. Additional experi- 
ments are necessary to confirm the response of the companion mode and to examine the 
nonlinear behavior for larger values of £. Measurements of the experimental vibration 
shape would also be worthwhile; the mode- shape data available to date are in substantial 
agreement with the radial deflection used in the analysis. 

It should be noted that the present theory imposes no restraints on the axial in- 
plane displacements at the ends of the cylinder. Imposing the boundary condition that 
the axial displacement u = 0 at the cylinder ends might strongly alter the nonlinear 
behavior, especially for high values of £ (that is, for "short" cylinders). In addition, 
it will be recalled that the present analysis loses accuracy for low values of the circum- 
ferential mode number n (for example, n < 5). This loss of accuracy is not expected 
to alter the qualitative behavior of the solutions, however. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 2, 1967, 

124-08-05-08-23. 
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APPENDIX A 


DISCUSSION OF THE CONTINUITY REQUIREMENT 
FOR THE CIRCUMFERENTIAL DISPLACEMENT 


Since the cylindrical shell is assumed to be complete and circular, the displace- 
ments, slope, moments, shears, and stresses must satisfy continuity requirements as 
noted in equation (3). In particular, the continuity requirement on the circumferential 
displacement v results in conditions which the functions w and F must satisfy. 

The purpose of this appendix is to derive and discuss these conditions for w and F. 

The continuity requirement for v is given by 

v(x,y,t) = v(x,y+27rR,t) (Al) 


This condition can be rewritten in the form 


py+2irR ~ 

v(x,y+27rR,t) - v(x,y,t) ■ 0 = J dy 


(A2) 


If an expression for 0v/8y is to be found in terms of w and F, it is necessary 
to examine the force-displacement relations. These relations are given by (see ref. 2) 

~\ 






9 2 F 




9x 2\Sxl 


8y w . 1 /9w V 
8y ” R + 2\8y/ 


= N,r = 


Eh 


!|V_W + l/9w\ +I/ 


8x2 y (l - U 2 ) fy R 2^/ 


&U l/8wV 
9x 2\8x/ 


(A3) 


J 


Eliminating 8u/ 9x from these egressions and then solving for 9v/8y gives 


8y _ 1 / 8 2 F 8 2 f\ w 1 / 8w\^ 

9y = Eh\Bx 2 ’ V &y2/ R ' 


(A4) 


Substituting this expression for 9v/ 8y into the continuity condition given by equa- 
tion (A2) results in the following relation between w and F: 



dy = 0 


(A 5) 


Equation (A5) represents a condition on w and F which must be satisfied if the 
circumferential displacement v is to be continuous. 
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When the expressions for w and F given in references 1 to 3 are substituted 
into equation (A5), the integral does not vanish. In other words, the functions for w 
and F given in these references will not satisfy the continuity requirement for the 
circumferential displacement v. 

In order to satisfy the continuity requirement, it is necessary to modify w or F, 
or both. For this reason, the deflection in the present report is chosen as given by equa- 
tion (7). When w and F from equations (7) and (8) are substituted into equation (A5), 
the integral vanishes as required. Thus, the continuity condition for v is satisfied by 
the w and F given in the present report. 
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APPENDIX B 


STABILITY ANALYSIS 


The purpose of this appendix is to discuss the details of the stability calculations 
and to derive the stability boundaries previously discussed. 


Stability of the Single-Mode Response 

When only the driven mode responds, the response is given approximately by 
(eq. (16)) _ . 

£ c (t) = A cos fir I 

e s (r) = o 


As a test of the stability of the response, small-perturbation terms 4> and tj 
are added to the approximate solution (Bl); that is, 

£ c (t) = A cos Or + $(r) | 

? s (t) = 0 + tj(t) [ 


Equations (B2) are substituted into equations (9) for j? and j? . The approximate 

c s 

solution given by equations (Bl) is assumed to satisfy equations (9) exactly, and the equa- 
tions governing (f> and r? are therefore homogeneous. Only first-order terms in the 
perturbations are retained, and the equations for <j b and 77 become 


1 + A 2 cos 2 Ot)^-^ - S2A^sin 2S2 r ^ + 1 - 
k 8 / dT 2 dT 


. o.a 2 cin 90. A 1 _ 3 C A 2 _ 3eyA 2 A 15 9^4 


^- + T^ A 


9e 2 6A 2 3eyA 2 5 2*T 4 
' 16 + 2 “ 2 C / 


- | e 2 6 A^)cos 2 Or + ^ cos 4 Sir 4> = 0 


-2 9 -4 

eyA 3e 2 6A 

2 8 


2-2 -2 9 — 4\ 9 -4 

3eO A .eyA e 2 6 A | , e 2 6 A ^ ^ A 

g + — ^ 2 — 7 C0S + — g — cos 4S2r 77 = 0 
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Equation (B3) is then simplified by the following transformations: 
Let 6 = S2t, so that 

d d dfl _ o d 
dr d0 dT d 6 

3eA 2 

Define Mj = “ — and substitute for <f>: 


8 


n6 fJ.i sin 26 dd 

4> = U(0)exp -j-± o - V 

^0 2(1 + iii cos* 6) 


Then, equation (B3) yields 


n 2 ^i + iii cos 2 #)^-^ + 


1 - pjO 2 cos 2 0 - 3eyA 2 cos^0 + 5«^6A^cos^0 


3p^f^sin 2 20 
4^1 + juj cos 2 e) 


U = 0 


(B5) 


This equation is further simplified by dividing by (l + fij cos^fl) and expanding the 
denominator for small values of /ip Retaining only first-order terms ^namely, 
of order eA 2 ) gives 


terms 


+ (M + 16K cos 20)U = 0 

dfl 2 


(B6) 


where 



(B7) 


16K = - 


Equation (B6) is now in the form of Mathieu's equation, which is discussed in 
detail by Stoker (ref. 12). The results of reference 12 are directly applicable to this 
analysis and can be used to determine approximate stability boundaries for the 
perturbation <p(r). 

The result of these calculations indicates that the perturbation <p is unstable 
within the region given by 
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<*»> 

iii 

and that smaller unstable regions exist near O = -, ~, -, . . . . 

it 3 4 

Equation (B4) for the perturbation 77 ( t ) can be treated in a similar manner if 
terms of order e 2 A 4 are neglected. Again, when the results for Mathieu’s equation 
are applied, the perturbation 77 ( 7 ) is found to be unstable within the region 

< b9) 

and in narrow regions near £2 = ^, . . . . 

Within the region denoted by the inequality (B 8 ), perturbations of the amplitude of 
the driven mode are unstable and jumps occur in the steady-state response. This region 
is the area where this analysis indicates that the jump phenomena of nonlinear vibrations 
occur. The region specified by the inequality (B9) represents an area where perturba- 
tions of the companion mode are unstable; in this region it is necessary to investigate 
coupled vibrations where both modes respond. For the case of coupled- mode vibrations, 
see equations (19) and (20) and the related discussion. The stability of the coupled-mode 
response is discussed in the following section. 


Stability of the Coupled-Mode Response 


When both modes respond, their amplitudes are given (eqs. (19)) approximately by 





(BIO) 


As a test of the stability of the response, the small perturbations <p and 77 are 
again introduced: 

“N 

£ c (t) = A cos S2 t + 0(t) 

> (Bll) 

S s (t) = B sin S2r + 77 ( 7 ) j 

These expressions are then substituted into equations (9) for £ c and £ g . The approxi- 
mate solution given by equations (BIO) is assumed to satisfy exactly equations (9), and 
only first-order terms in the perturbations are retained. This procedure results in two 
coupled stability equations for 0 and 77 , namely, 
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^1 + cos 2 ^ sin 20 + (p + jjq j O 2 ^sin 2 6 ~ 2 cos 2 0^ + [i^l^cos 29 

(3A 2 cos 2 0 + B 2 sin 2 0)+ e 2 6(5A 4 cos^0 + 6A 2 B 2 sin 2 0 cos 2 0 + B^cos^jJ (f> 


- ey 


and 


2 r 

+ AB sin 6 cos 9 ^-5 + ^ fJAB cos 2 6^p- - [^(s2 2 AB sin 6 cos 0 ) 

8 dr 2 4 dr [8 ' ' 

+ eyAB sin 20 - 2e 2 6AB sin 20 (a 2 cos 2 0 + B 2 sin 2 0^| tj = 0 

^1+ ju-2 cos ^^)^2 + ^2^ s * n dr + ^ + |^i^ cos 20 + M 2^( cos ^® " 2 sin 2 0) 

- ey(A 2 cos 2 0 + 3B 2 sin 2 0) + e 2 e(A^cos^0 + 6A 2 B 2 sin 2 0 cos 2 0 + 5B 4 sin^0)] r\ 

+ 4r AB sin 0 cos 0 S2AB sin 2 0 ^ + P^- J2 2 AB sin 0 cos 0 

8 dr 2 ^ dr L8 


(B12) 


+ eyAB sin 20 - 2e 2 6AB sin 20 (a 2 cos 0 + B^sin 2 ©) <f> = 0 


(B13) 


where the substitutions 


0 = S2t 


Mj = 

M2 = 


3eA 2 


8 

3eB 2 

8 


have been used to simplify the results. 

Equations (B12) and (B13) are coupled equations with periodic coefficients, and 
exact solutions to them are not known. However, approximate results for the stability 
boundaries of these equations can be obtained by using the method of slowly varying 
parameters (ref. 19). 
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In this method, the perturbations <fi and r\ are assumed in the form 

<£(t) = 4>^(t ) cos SIt + (^(Tjsin fi-rj 

t] (r) = ?Jj(t)cos Sir + r72(T)sin Sir 


(B14) 


where <^2> Vp and 772 are assumed to be slowly varying functions of r. Then, 

the derivatives d<p/dr, d^/dr 2 , d rj/dr, and d^rj/dr^ are replaced by 

^ sin Or + cos ^ T 


j 2 . 2 o d0< d<^>2 

2— 2 = -fi cos Or - 0^02 sin Ot - O — pj=- sin Or + O cos S2r 


dr 2 

S~°»1 sin J2 t + O772 cos fir 

j|2 n 2 dlj - | d7Jn 

= -fo?! cos S2r - O T72 sin Or - f2 sin Or + SI cos J2r 


> (B15) 


together with the auxiliary conditions 

d <t> 1 


's 


dcp 

, cos Sir + — sin Ot = 0 

dr dT 


drjj dr? 2 

cos J2r + sin J2 t = 0 

dr dT 


(B16) 


These expressions for the derivatives are then substituted into equations (B12) and (B13), 
and the auxiliary conditions (B16) are used. The resulting equations are multiplied by 
cos J2r and averaged by integrating t from 0 to 2tt. In the integration, the variables 
$1 through ri 2 are replaced by their average values and denoted by a bar - i.e., 

<pp . . . 77 2 - This procedure yields two equations for 0 2 , V-^, and ij 2 * Two 

more equations are obtained by returning to equations (B12) and (B13), multiplying by 
sin f It, and performing a similar integration. The final result is four linear differential 
equations which can be put in matrix form as follows: 

[?]<? = Wf? (B17) 
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where 


r~~\ 

4>1 

02 



^2 


and the matrices [p] and [T] contain elements that are constants which depend upon 
the steady- state solution through the values of A, B, O, e, y, and 6. 

Equation (B17) is solved by using cp = <PqG^ T and leads to an eigenvalue problem 
for the X: 

|P-XT|=0 (B18) 

The matrices [Pj and [T] are nonsymmetric, and equation (B18) yields complex 
roots. If any of the roots have a positive real part, then the corresponding perturbations 
increase exponentially with time. In this case, the associated steady-state solution is 
said to be unstable; conversely, if none of the roots of equation (B18) has a positive real 
part, the solution is said to be stable. 

Stability of the steady-state solution given by equation (BIO) was investigated by 
substituting representative values of A, B, and into equation (B18) along with the 
associated parameters e, y, and 6. For each case, the roots were examined to 
determine whether or not they had a positive real part. In this manner, the stability of 
the solution plotted in figure 5 was determined; the results of the calculations are dis- 
cussed in the main text. 
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